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Abstract 

We study large N SU (N) Yang- Mills theory in three and four 
dimensions using a one-parameter family of supergravity models which 
originate from non-extremal rotating D-branes. We show explicitly 
that varying this "angular momentum" parameter decouples the Ka- 
luza-Klein modes associated with the compact D-brane coordinate, 
while the mass ratios for ordinary glueballs are quite stable against this 
variation, and are in good agreement with the latest lattice results. We 
also compute the topological susceptibility and the gluon condensate 
as a function of the "angular momentum" parameter. 
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1 Introduction 



Generalizing the conjectured duality |IJ between large N superconformal field 
theories and superstring or M theory on anti-de Sitter (AdS) backgrounds, 
Witten proposed an approach to studying large N non-supersymmetric the- 
ories such as pure QCD using a dual supergravity (string theory) description 
[0]. The basic idea is to start with d + 1 dimensional superconformal field 
theories at finite temperature - thus breaking the superconformal invari- 
ance - and obtain a d dimensional non-supersymmetric gauge theory at zero 
temperature by dimensional reduction in the Euclidean time direction. The 
AdS space is then replaced by a certain limit of the Schwarzschild geometry 
describing a black hole in an AdS space. 

When the curvature of the space is small compared to the string scale (or, 
in the case of M theory, Planck scale), supergravity provides an adequate ef- 
fective description that exhibits a qualitative agreement with pure QCD in 
three and four dimensions 0, [|. The supergravity limit of string theory (i.e. 
infinite string tension, a' — > limit) corresponds to the strong coupling limit 
of the gauge theory (A = gy M N 3> 1), with 1/A playing the role of a'. In the 
approach of 0], the gauge theory has an ultraviolet cutoff proportional to the 
temperature T; the supergravity approximation should describe the large TV 
gauge theory in the strong coupling regime with a finite ultraviolet cutoff. 
This is analogous to a strong coupling lattice gauge theory with lattice spac- 
ing ~ 1/T |. In the limit that the ultraviolet cutoff is sent to infinity, one 
has to study the theory at small A, and the supergravity description breaks 
down. To calculate the spectrum in this regime, a better understanding of 
string theory with Ramond-Ramond (R-R) background fields is required. 

Glueball masses in the supergravity approximation have been computed 
in |4], [|. The Witten model contains in addition to the glueballs certain 
Kaluza-Klein (KK) particles with masses of the order of the glueball masses. 
These KK modes do not correspond to any states in the Yang-Mills theory, 
and therefore they should decouple in the "continuum" limit. The KK states 
do not decouple with the inclusion of the leading 1/A corrections ||. Al- 
though such states can decouple in a full string theory treatment, there may 
be generalizations of the Witten model that have a more direct connection 
with the continuum gauge theory already at the supergravity level (at least 
in the sector of states with spin < 2 that can appear in a supergravity de- 
scription) . A similar situation arises in lattice gauge theory. It is well known 
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that the action one starts with has a significant effect on the speed at which 
one gets to the continuum limit. One can add to the lattice action defor- 
mations which are irrelevant in the continuum and arrive at an appropriate 
effective description of the continuum theory while having a larger lattice 
spacing (such a deformed action is called an "improved" lattice action). A 
similar strategy in the dual supergravity picture would correspond to a suit- 
able modification of the background metric, so as to have an appropriate 
effective description of the gauge theory while still having a finite ultraviolet 
cutoff. An important test of the proposal is to check that the KK modes 
in the supergravity description that do not correspond to gauge degrees of 
freedom are heavy and decouple, and at the same time the infrared physics is 
not significantly altered. In this paper we make the first step in this direction 
by examining a generalization of the Witten model that has an additional 
parameter. 

A more general approach to the conjectured correspondence between 
gauge theories and M-theory requires the investigation of supergravity com- 
pactifications which asymptotically approach anti-de Sitter backgrounds, e.g. 
AdS$ x S 5 or AdS-j x S 4 (see e.g. |7], ||). There exist a few supergravity back- 
grounds that generalize the Witten model and which are regular everywhere. 
These are essentially obtained by starting with the rotating version of the 
non-extremal D4 brane background (or rotating D3 brane background, in 
the case of QCD 3 ) and taking a field theory limit as in [p]]. These models 
were investigated in |§. The deformation of the background proposed in || 
is parameterized by an "angular momentum" parameter (the supergravity 
background is actually static, with the Euclidean time playing the role of an 
internal angle). In this paper we determine numerically the scalar and pseu- 
doscalar spectrum of these models as a function of the angular momentum 
parameter and compare the results to those obtained by lattice calculations. 
We also compute the gluon condensate and the topological susceptibility 

The paper is organized as follows: Section 2 is devoted to the study of 
supergravity models for pure QCD in 3 + 1 dimensions. The models are 
described in Section 2.1. In Section 2.2 we compute the scalar glueball mass 
spectrum and analyze its dependence on the angular momentum. In Section 
2.3 we calculate the mass spectrum of some KK modes. It will be shown that 
the KK modes associated with the compact D-brane coordinate decouple 
as the angular momentum parameter is increased. This, however, is not 
the case for the SO(3) non-singlet KK modes with vanishing U(l) charge 



2 



in the compact D-brane coordinate. In Section 2.4 we compute the gluon 
condensate from the free energy associated with the supergravity background. 
In Section 2.5 we compute the topological susceptibility and its dependence 
on the angular momentum parameter. Section 3 contains an analogous study 
for the case of QCD in 2 + 1 dimensions. The conclusions are similar in both 
cases, and they are summarized in Section 4. 



2 QCD in 3 + 1 dimensions 

2.1 Supergravity Models for QCD 4 

One way to construct non-supersymmetric models of QCD based on super- 
gravity is to start from the non-extremal D4 brane metric, and view the 
Euclidean time coordinate as an internal coordinate compactified on a circle 
of radius {2-kTh)^ 1 [0. Possible generalizations of this proposal are con- 
strained by the no-hair theorem, which implies that the most general regular 
manifold with only D4 brane charges (and an isometry group containing R 4 ) 
is given by the rotating version of the non-extremal D4 brane, which has two 
additional parameters representing angular momenta in two different planes. 
The Euclidean version of this metric (related to the rotating M5 brane metric 
by dimensional reduction) was used in || to construct models for QCD with 
extra parameters. Here we investigate in detail the case when there is one 
non- vanishing angular momentum, parametrized by a. The field theory limit 
of the Euclidean rotating M5 brane with angular momentum component in 
one plane is given by the metric || 



+ (vriV) 2 / 3 (dO 2 + ^ sin 2 Odip 2 + i cos 2 9dQ 2 2 

2a 2 f/«? . 2 „, , A 4(ttA0 2 /W 
sm 2 6drdip + 



Z7 4 A(ttJV) 1/2 / u 2 (l-^-% 



(2.1) 



where xi, . . . ,x§ are the coordinates along the M5 brane where the gauge the- 
ory lives, U is the "radial" coordinate of the AdS space, while the remaining 
four coordinates parameterize the angular variables of S 4 , and where we have 
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introduced 

a 4 cos 2 9 ~ a 4 
A = 1 — , A = 1 . (2.2) 

Dimensional reduction along x$ (which will play the role of the "eleventh" di- 
mension) gives N rotating non-extremal D4 branes, which in the low energy 
regime should be described by SU(N) Yang-Mills theory at finite tempera- 
ture T H , perturbed by some operator associated with the rotation. The 3 + 1 
dimensional SU (N) Yang-Mills theory at zero-temperature can be described 
by making x\ — > — ixo, and viewing r as parameterizing a space-like circle 
with radius Rq — (27rT^) _1 , where fermions obey anti-periodic boundary 
conditions. At energies much lower than 1/Ro, the theory is effectively 3 + 1 
dimensional. Because of the boundary conditions, fermions and scalar parti- 
cles get masses proportional to the inverse radius, so that, as -R - ► 0, they 
should decouple from the low-energy physics, leaving pure Yang-Mills theory 
as low-energy theory. 

The gauge coupling g\ in the 3 + 1 dimensional Yang-Mills theory is given 
by the ratio between the periods of the eleven-dimensional coordinates x 5 and 
r times 2n. It is convenient to introduce ordinary angular coordinates 9±, 
and 9 2 which are 27r-periodic by 

r = R 9 2 , x 5 = ^-R 9 1 = ^-R 9 1 , (2.3) 
R = (27rTf/) 1 = A , A = 4 U \ 4 , (2.4) 

Silo Ujj — gGT 

where uh is the location of the horizon, and we have introduced the 't Hooft 
coupling 

the coordinate u by U — 2(tcN) 1 ^ 2 u, and rescaled a — > 2(nN) l l 2 a. By 
dimensional reduction in 9±, one obtains the metric 



<4a = 

3u 



Au ( — dx + dx 1 + dx 2 + dx 3 ) 



AA 2 2 u§ 4 ^ 



o - 



A 1 An 2 An 2 

+ d ^ £ gin 2 ^ 2 1 cog 2 ^ Q 2 _ 40^o ^2 m d 

A A 2 3w 4 A 



,(2.6) 
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with a dilaton background 



6 ~ 27 ujj ^2 • ^-0 

With this normalization, the metric reduces to Eq. (4.8) of ref. after 
setting a = 0. The string coupling is of order 1/N, and the metric has 
become independent of N, which is consistent with the expectation that in 
the large N limit the string spectrum should be independent of N. The 
metric is regular, and the location of the horizon is given by 

u% - a 4 u 2 H - ul = , (2.8) 

i.e. 



u 



2 fl4 




+ M, 7 = 3 1 + - • (2.9) 



Note that for large a, one has the approximate expressions 

(u% is always real). This shows, in particular, that the radius Rq = A/(3uq) 
can be made very small by increasing a/u . This is essentially the mechanism 
that will make the corresponding KK modes decouple at large a/w . At small 
a, such KK states have masses of the same order as the masses of the lightest 
glueball states. 

The string tension is given by H] 



4, . o _uS 
4 



String excitations should have masses of order a 1 ! 2 . The spin < 2 glueballs 
that remain in the supergravity approximation - whose masses are deter- 
mined from the Laplace equation - have masses which are independent of A. 
The supergravity approximation is valid for XA ^> 1 so that all curvature 
invariants are small 0. In this limit the spin > 2 glueballs corresponding to 
string excitations will be much heavier than the supergravity glueballs. 
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2.2 Spectrum of glueball masses 

The glueball masses are obtained by computing correlation functions of gauge 
invariant local operators or the Wilson loops, and looking for particle poles. 
Following |7|, gj, correlation functions of local operators O are related at large 
N and large <7y M -/V to tree level amplitudes of supergravity. The generating 
functional for the correlation function of O is the string partition function 
evaluated with specified boundary values <^o of the string fields. When the 
supergravity description is applicable we have 

/ d 4 x<p (x)O(x)^ _ e -Isa(vo)^ (2.12) 

where Isg is the supergravity action. 

The spectrum of the scalar glueballQ ++ is obtained by finding the nor- 
malizable solution to the supergravity equation for the dilaton mode $ that 
couples to TrF 2 , which is the lowest dimension operator with ++ quantum 
numbers. 

The equation for $ reads 

d^e^g^d^} = 0, (2.13) 

where is the string frame metric. 

We look for ^-independent solutions of the form $ = x{ u ) elkx - One 
obtains the equation 

- a\ 2 - uDx'M] = -M 2 X (u) , M 2 = -k 2 , (2.14) 

u 6 

where the eigenvalues M are the glueball masses. The solution of this (ordi- 
nary) differential equation has to be normalizable and regular both at u — > uh 
and u — > oo. The eigenvalues of this equation can be easily obtained numer- 
ically M by using the "shooting" method. One first finds the asymptotic 
behavior of xi u ) f° r large u, and then numerically integrates this solution 
back to the horizon. The solutions regular at the horizon will have a finite 
derivative at uh- This condition will determine the possible values of the 
glueball masses M. The results of this numerical procedure are presented in 



Table [2.1| and in Figure EO. One can see from Figure |2.1| that the ratios of 



2 In the following we will use the notation J PC for the glueballs, where J is the glueball 
spin, and P, C refer to the parity and charge conjugation quantum numbers respectively. 
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the masses of the excited glueball states compared to the ground state are 
very stable with respect to the variation in the parameter a, even though 
both quantities themselves grow like M 2 ex a 2 for large a. The asymptotic 
value of the mass ratios is taken on very quickly, a/uo ~ 2 is sufficiently 
large to be in the asymptotic region. A priori one could have expected that 
these mass ratios may change significantly when a is varied. This leads one 
to suspect that there is a dynamical reason for the stability of the ratios of 



masses. 



state 


lattice, N = 3 


supergravity a = 


supergravity a — > oo 


0++ 


1.61 ±0.15 


1.61 (input) 


1.61 (input) 


o++* 


2.8 


2.55 


2.56 


Q++** 




3.46 


3.48 


qH — h*** 




4.36 


4.40 



Table 2.1: Masses of the first few ++ glueballs in QCD 4 , in GeV, from 
supergravity compared to the available lattice results. The first column 
gives the lattice result [HJ PI , the second the supergravity result for 
a = while the third the supergravity result in the a — > oo limit. The 
rrj do note quote an error on the preliminary lattice 

oo in the 



authors of ref. 
result for ++ 



Note that the change from a 
supergravity predictions is tiny. 



to a 



Let us now consider the _+ glueballs. The lowest dimension operator 
with _+ quantum numbers is TrFF. On the D4 brane worldvolume, the 
field that couples to this operator is the R-R 1-form A^. In order to find the 
0~ + glueball masses we have to solve its equation of motion 



d u [^ggWg™(d p A a -d a A p ) 
Consider solutions of the form 

Ae 2 = Xe 2 {u) e" 







1,...,10 



Ak-x 



A^ = if fi^9 2 



Plugging this into ( 2.16|) , we obtain 

which reads 
1 , , 



(2.15) 
(2.16) 
(2.17) 



aV - u 6 )d u [u%u 4 - a 4 ) X ' e , 2 (u)} = ~ M V - a 4 ) X e 2 (u). (2.18) 
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Figure 2.1: The dependence of the ratio r = of the masses of the 

first excited (0 ++ *) glueball state to the lowest ++ glueball state on 
the parameter a (in units where uq = 1). The ratio changes very little 
and takes on its asymptotic value quickly. 



For a = 0, it yields Eq. (2.9) of ||12|| , as required. When a ^ there are no 
solutions of the form ( |2.16| ). The reason is that the gg 2lfi component of the 
metric ( |2.6| ) is non- vanishing for a ^ and, as a result, the tp component 
of the Maxwell equation is not satisfied automatically (solutions contain a 
non- vanishing component A^). 

We will work in the approximation a/uo ^> 1. In this approximation the 
non-diagonal gg 2V part of the metric can be neglected, and there are solutions 
of the form ( |2.16|) . Effectively, we can solve (|2.18|) in the limit a ^> Uq. We 
must however keep in mind that we need Uq ^ to regularize the horizon, 
and the actual limit that is taken is a/uo large at fixed Uq (so that curvature 
invariants are bounded from above and they are small for sufficiently large 
't Hooft coupling A). 

The mass spectrum from ( |2.18p can be obtained using a similar numerical 
method as for the ++ glueballs. The dependence of the lightest 0~ + glueball 
mass on a is presented in Figure |2j|, whereas the h glueball mass spectrum 
in Table [2.2| . Note that while masses ratios are fairly stable against the 
variation of a (they again grow like M 2 oc a 2 ), the actual values of the mass 
ratios compared to ++ increase by a sizeable (~ 25%) value. The change 
is in the right direction as suggested by recent improved lattice simulations 
13] . The mass of the second 0~ + state also increases and is in agreement 
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with the new lattice results fT3|| . 



state 


lattice, N = 3 


supergravity a = 


supergravity a — > oo 


o-+ 


2.59 ±0.13 


2.00 


2.56 


-+* 


3.64 ±0.18 


2.98 


3.49 


Q 




3.91 


4.40 


Q — h*** 




4.83 


5.30 



Table 2.2: Masses of the first few 0~ + glueballs in QCD 4 , in GeV, 
from supergravity compared to the available lattice results. The first 
column gives the lattice result, the second the supergravity result for 
a = while the third the supergravity result in the a — > oo limit. Note 
that the change from a = to a = oo in the supergravity predictions 
is sizeable, of the order ~ 25%. 



We can directly compare the ratio of masses of the lowest glueball states 
h and ++ with lattice results [f[0|, O, [03]. Since one of the largest errors 



in the lattice calculation of glueball masses comes from setting the overall 
scale?] the ratios of masses are even more accurately known from the lattice 
than the masses themselves. Using the lattice results UTT , 13, [14[| in the more 
accurate "lattice units" r : 



r M 0+ + = 4.33 ±0.05, r M -+ 
we find: 



6.33±0.07, r M -+* = 8.9±0.1, (2.19) 



\a=0 



0++ * supergravity 



'hi 



' M, 



o- 



\ a— >oo 



0++ / 



supergravity 



lattice 



1.24 
1.59 

1.46 ±0.03 



(2.20) 



Mr, 



a=0 



0++ ' supergravity 



0++ * supergravity 



J 0++ / lattice 



1.85 
2.17 

2.06 ±0.05 



(2.21) 



3 We thank M. Peardon for emphasizing this point to us. 



9 



1.6 
1.55 



/ 

1.5 [ / 



1.45 
1.4 
1.35 

1.25 



i 



I 

I 

I 

I 

a 
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Figure 2.2: The dependence of the ratio r = on ^ ne parameter a 

(in units where Uq = 1). The change in the ratio is stable against the 
variation of a, however it increases by about 25% while going to a = oo. 
The change is in agreement with lattice simulations. As explained in 
the text, this figure is reliable only for the regions a < «o or a > ii 
which are shown in the plot with a solid line, while for the intermediate 
region denoted by a dashed line there are corrections due to the non- 
vanishing off-diagonal component of the metric go 2ip . 
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One can see that taking the a — > oo improves the agreement between the 
supergravity and lattice predictions significantly. One should however keep 
in mind that the supergravity results presented here are for the limit N — > oo 
and A — > oo, while the lattice results are for N = 3 and A small. [] Direct 
lattice calculations for the large N limit are have just started to become 



available [15], however no reliable direct estimate for the mass of the + is 
known yet. 

2.3 Masses of Kaluza-Klein states 

In the supergravity approximation the a = model contains additional light 
KK modes in the spectrum whose masses are of the same order as those of 
the glueball states f6j. In this section we investigate whether the additional 
parameter of the model considered here can be tuned to decouple the KK 
modes already at the supergravity tree-level. In the following, it will be 
shown that this is indeed the case for the KK modes wrapped around the 82 
direction, which become very heavy for a/u ^> 1. We thus look for solutions 
of the dilaton equation (|2.14|) of the form 

$ = x (u) e ik - x e ind2 . (2.22) 

One finds the following equation: 

±d u Uu°-aW-4) X '(u)) = (~M 2 + ^ / 4 ~/ M a ) X(u) , (2.23) 
u \ j \ -rto w 4 — a 4 % J 

where Rq is given in Eq. Q2.4|) . This generalizes fl2.14j ) to the case n^O. We 
want to compare M = M(n = 0) with Mkk = M(n = 1). The question is 
how Mq/Mkk behaves as a function of a (we can set Uq = 1). The extra term 
proportional to l/Rl gives a positive contribution to the mass, so that M KK 
should increase as M\ K oc a 8 as a is increased (the KK radius Ro = A/3uo 
shrinks to zero as a —>■ 00). Thus one expects that Mo/Mkk oc 1/a 3 — > 
as a increases. The numerical values of the masses of these KK modes are 



displayed in Table |2.3| and Figure 2.3 . Note that the numerical evaluation of 



the masses of these KK modes becomes more and more difficult as a increases. 



4 For example, in Ref. (jlll the results are extrapolated to A = from calculations in 
the region 7.5 < g 2 N < 10. 
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This is because the term with 1/Rq causes an overall shift of the masses, while 
the splittings between the excited KK modes still remain of the same order as 
for the ordinary glueballs. As a result, the solutions become more and more 
quickly oscillating as a increases, making numerical treatments increasingly 
difficult. For this reason we display only values up to a/u = 3. 



state 


value for a = 


value for a = 3 


KK 


2.24 


20.25 


KK* 


3.12 


20.37 


KK** 


4.01 


20.52 




4.89 


20.72 



Table 2.3: Masses of the KK modes which wrap around the 62 circle 
and have no corresponding states in QCD 4 , in GeV. The first column 
gives the masses for a = while the second the masses for a = 3. Note 
that even for a = 3 these states are heavier by a factor of 10 than the 
++ glueball mass and are effectively decoupled from the spectrum even 
in the supergravity limit. 

Above we have demonstrated that the KK modes which correspond to 
states that wrap the 62 direction are effectively decoupled from the spectrum 
even in the supergravity approximation. However, there are other KK modes 
in this theory, and one would like to know whether these are decoupled as 
well. The reason for the decoupling of the modes on 9 2 is clear: the radius of 
this direction shrinks to zero when a — > 00. However, the radii of the other 
compact directions do not behave similarly. Therefore it is reasonable to ex- 
pect that these states will not decouple at the level of supergravity from the 
spectrum (but they could decouple once string theory corrections are incor- 
porated). We now demonstrate by explicit calculation of the corresponding 
mass spectrum that this is indeed the case. 

Consider non-singlet modes which are independent of 82, of the form 
f(u)e lk ' x cos if sin 9. This corresponds to spherical harmonics on S A with an- 
gular momentum 1 = 1. Plugging this ansatz into the dilaton equation ( |2.13j ) 
we find that f{u) satisfies the equation (uq = 1) 

u d u [(u — a u — u)f (u)] — f(u) I k + \— — =0. 2.24 



12 




Figure 2.3: The dependence of the ratio r = M 0+ * of the the lowest 
++ glueball state compared to the KK mode wrapping the #2 circle on 
the parameter a in units where uq = 1. This KK mode decouples very 
quickly from the spectrum even in the supergravity approximation. 



The results of the numerical analysis of the eigenvalues are presented in 



Table £]4] and Figure |2]J. One can see that these states do not decouple 
from the spectrum at the supergravity level, instead their masses remain 
comparable to the ordinary glueball masses. 



state 


value for a = 


value for a — > 00 


KK 


2.30 


2.84 


KK* 


3.29 


3.80 


KK** 


4.23 


4.74 




5.15 


5.65 



Table 2.4: Masses of the KK modes corresponding to I = 1 angular 
momentum on the S 4 , in GeV. The first column gives the masses for 
a = while the second the masses in the a —>■ 00 limit. Note that the 
change from a = to a = 00 in the supergravity predictions is not 
sufficiently large in order to decouple these particular states from the 
spectrum in the supergravity limit. 
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r 




a 



Figure 2.4: The dependence of the ratio r = M K * of KK modes (cor- 
responding to spherical harmonics with I = 1 on S 14 ) compared to the 
lowest ++ glueball state on the parameter a in units where uq = 1. 
This KK mode does not decouple from the spectrum in the supergravity 
approximation even in the a —>■ oo limit. 



2.4 Free energy and Gluon condensation 

The standard relation between the thermal partition function and free energy 
Z(T) = exp(— F/T) relates the free energy associated with the supergravity 
background to the expectation value of the operator TrF^ u . This relation 



was exploited in |12| to obtain a prediction for the gluon condensate in the 
Witten (a = 0) supergravity model. Let us now derive the corresponding 
supergravity result for general a. From the rotating M5 brane metric (given 
in Eq. (3.1) of ||), one can obtain the following formulas for the ADM mass, 
entropy and angular momentum (see also WQ, pi): 



Ma.dm 
S 

Jh 



W 9 n 3 . 2 
— — — — 2m(l + - smh a) 
AnG N 4 



v(n 4 



8vr 2 



v 5 v(n 4 

V 5 V(tti 
4ttG n 



-Imru cosh a , 
-ml cosh a , 



(2.25) 
(2.26) 
(2.27) 



G 



K 



N 
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8tt 



p > 



(2.28) 
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where Gat is Newton's constant in 11 dimensions, and lp is the 11 dimensional 
Planck length. The (magnetic) charge N is related to a and m by 

2m cosh a sinh a = nNlp . (2.29) 

The Hawking temperature and angular velocity are given by 

3r?r + I 2 lr H , 

T H = H , Q H = ~ ^r- ■ 2.30 

Snm cosh a 2m cosh a 

These quantities satisfy the first law of black hole thermodynamics: 

dM ADM = T H dS + n H dJ H . (2.31) 

We are interested in the field theory limit lp — > 0, obtained by rescaling 
variables as follows 

r = u%{A-kN) , m = ^u 6 l 9 P (^N) 3 , I = ia%{AixN) . (2.32) 



We get 



E = M ADM - M extremal = — V 5 N\ 6 , (2.33) 



Th = ^, S=^N 3 u 2 h u 3 , (2.34) 



Q H = i 2 -^- , J H = i^V,N 3 a 2 u 3 , (2.35) 

Uq 07T 



with 



A - , u% - a* = 4 • (2.36) 

6u H — a 4 u H 



The free energy is given by 



E-T H S- n H J H = ~^N 3 u 6 . (2.37) 

07T° 



For a = this reproduces the result of |18[ («o = 2ttATh/S). The M5 brane 
coordinate 25 is compactified on a circle with radius RoX/N, given by Eq. 
so that 

^5 = 7^7 • (2.38) 
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The gluon condensate is then given by 

4#ym Vi± H LtTX 

For a = this reduces to the corresponding result in [T^] (setting A = 1 and 
u = 2ttT h /3). Expressing u in terms of the string tension ( |2.11|) we obtain 

Note that this relation is independent of a (in particular, it applies to the 
a = case as well). It has the expected dependence on N, and a simple 
dependence on A. 

2.5 Topological Susceptibility 

The topological susceptibility Xt is defined by 

Xt = J d 4 x (Tr FF (x)Tr FF(0)) . (2.41) 

The topological susceptibility measures the fluctuations of the topological 
charge of the vacuum. At large N the Witten-Veneziano formula |19], re- 
lates the mass m v r in SU (N) with Nf quarks to the topological susceptibility 
of SU(N) without quarks: 

m l = -JTXt • (2-42) 

The effective low-energy four dimensional brane theory contains the cou- 
pling 

J d A xd6 2 Ae 2 Tr F, u F Xa e^ , (2.43) 

where Ag 2 is the component along the coordinate 9 2 of the R-R 1-form A^. We 
will consider zero mode (M 2 = 0) configurations where Ag 2 is independent 
of the world-volume coordinates. Comparing to the standard Yang-Mills 
coupling, 

' ' d 4 x 9 Tr FF , (2.44) 



1 



16vr 2 
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one obtains the relation 



2- 



d9 2 A d2 = 2uA e2 . (2.45) 

The action of the R-R 1-form is given by 

1 = T^T I d w x^-(d,A v - dvA^Av, - dv.A^gW' g w ' . (2.46) 

As discussed in Sect. 2.2, in the approximation that a/u is either very large 
or very small, the metric is diagonal and there are zero mode solutions of the 
form Ag 2 = Ag 2 (u), A^ = , /i ^ 9 2 . The action reduces to 



, dAe 2 (u) 2 uu e 2 e 2 
1 du ' y y 



(2.47) 



Using Eq. (|2.6|) and integrating over the angular coordinates, this becomes 



2 7 n e A 2 X 3 v 
27mo k io 



u 3 (u 4 — ct 4 )( 



du 



The equation of motion is then given by 

d u [u 3 (u*-a 4 )d u Ae 2 ] =0 . 

Whence 



<Me 2 = 6A™C(a 



tip. 



uHu A - a 4 ) ' 



Ae 2 = A 



1 + 6C(a 



u 2 — a 2 



(2.48) 



(2.49) 



(2.50) 



(2.51) 



V2a 4 w 2 4a 6 * u 2 + a 2 . 
The integration constant C (a) will be fixed by assuming that Ag 2 (u) vanishes 
at the horizon This gives 

1 



C(a) 



3m§ 3u% u 2 H + a z 
h& + 2^ g ^| 



(2.52) 



The other integration constant is related to the ^-parameter by (|2.45|) , 
27rA^ = 8. Note that in the limit a = one gets C(0) = 1 and v4 e2 (-u) = 

^(i-S)- 
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Using kIq = k^ 1 /27t = 2 6 7r 7 and uq = 2'kATh/3, we obtain 

I = 9*V^A*C(a)\^ H . (2.53) 

The topological susceptibility ( |2.41|) can then be obtained by differentiating 
twice with respect to 6: 

327T 

Xt = ^A*C(a)\*T* , (2.54) 
or, in terms of the string tension ( |2.11|) , 

* = ^ 2 ■ (2-55) 

The ^-dependence of the vacuum energy of the form 6 2 is the result antici- 
pated in |2lJ for the a = model, and Eq. (|2.53| ) shows that it holds for large 
a too. In the large N limit, this must be the case for consistency |21|. For 
a = 0, one has A = 1 = C(a), and Eq. (|2.54|) reproduces the result obtained 
in (12j . In the large a limit we have (see Eq. ( p,10| )) 

« n..6?l a , (2-56) 



9«ologf- 



so that 



"0 

I ,: S "1! 



Xt « 7^—^— = _- 5 _ r A^. (2.57) 



18tt 3 a 2 log 72tt 3 u§ log ^ 
This decreases if we increase a/uo at fixed A and uq. 



3 QCD in 2 + 1 dimensions 

3.1 Supergravity Models for QCD% 

Analogous models for QCD 3 can be obtained by starting with the Euclidean 
rotating D3 brane (Eq. (3.16) in || with x — > —ir, I — > il) and taking 
a! —>■ by rescaling variables as follows: 

r = Vol , 2m = Ufa' 4 , I = aa' . (3.1) 
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In the limit a' — > at fixed U, a, Uq we obtain 

ds\ lB = a'A^ 2 (h dT 2 + dx\ + dx\ + dx\) H ^— — 

1 - ■ ■ " ' En_ 



I' 4 

xa 

C/ 2 (7 4 - 



v^iV v 1 ' dy £/ 2 (l-^|- 

A . 2 „ , 2 cos 2 ,„ \ 2a?7 2 

- sm flay? H — 

•o A 



+ V / 4^ivU 2 + ^ sin 2 0dV 2 + ^-dtl 2 ) - sin 2 fldrdp 

V An An / U An 



(3.2) 



where 



a 2 



{7n » O? COS 2 r 

^0 = 1-^, Ao = l--^-, A ° = 1 -^' ( 3 ' 3 ) 

dfi 2 = d^i + sin2 + cos2 ^i#3 ■ ( 3 -4) 

The theory describes fermions with anti-periodic boundary conditions on 
the circle parameterized by r, which has radius (2ttTh)~ 1 with 

Th = , B = — j— -= — , (3.5) 

7r£? u H (2u z H — a 2 ) 



u 2 H = -{a 2 + ^a± + Aui) . (3.6) 

For convenience, we have rescaled variables by U — (Airg s N) l l 2 u , a — > 
(4vrc/ s iV) 1/2 a. At energies much lower than T# the theory should be effec- 
tively 2 + 1 dimensional (with xn = ix^ playing the role of time). The gauge 
coupling of the 2 + 1 dimensional field theory is given by 

5ym 3 = 9ym 4 t h , 9ym 4 = 2?r 9s ■ (3.7) 
A ss ^ = 9sN ^ (3.8) 

Z7T 7T.D 

In this model Wilson loops exhibit an area-law behavior with string tension 

1 

2^ 



a = ~ x /4ng s N u\ = \f\B uj 2 . (3.9) 



This can be obtained by minimizing the Nambu-Goto action P2"| , and it is 
essentially given by the coefficient of Yh=i dx 2 at the horizon times l/(2n) 0, 
|23|, ^4j. In the limit of large 't Hooft coupling, the light physical states are the 
supergravity modes, whose masses can be determined from the equations of 
motion of the string theory effective action. In the next sections we calculate 
the mass spectrum of the light physical states and of KK modes. 
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3.2 Spectrum of glueball masses 



In order to find the spectrum of the ++ glueball states one has to consider 
the supergravity equation for the dilaton mode $ that couples to the operator 
TrF 2 

dvy/ggrd v * = 0, (3.10) 

evaluated in the above background. For functions of the form $ = x{U) e lk ' x , 
we obtain 



tt [(u 4 - Wq - aV)ux'(u)] + M 2 u X {u) = 0, M = -k 



(3.11) 



The eigenvalues of this equation can again be determined numerically. The 

~~ J] gives the depen- 



results are presented in Table [3J] and Figure |3J]. Figure 
dence on a of the mass ratio of the first excited ++ glueball state compared 
to the ground state ++ . One obtains a very similar behavior to the case of 
QCD 4 , that is the mass ratio changes very little, and takes on its asymptotic 
value quickly. The comparison to the available lattice results [2S| are given 
in Table 13.11 



state 


lattice N — > oo 


value for a = 


value for a — > oo 


0++ 


4.065 ± 0.055 


4.07 (input) 


4.07 (input) 


++* 


6.18 ± 0.13 


7.03 


7.05 


Q++** 


7.99 ± 0.22 


9.93 


9.97 


qH — |-*** 




12.82 


12.87 



Table 3.1: Masses of the ++ glueballs and their excited states in QCD 3 . 
The first column gives the lattice results extrapolated to iV — > oo, the 
second column the supergravity results for a = and the third column 
the supergravity limit a — > oo. The lattice results are in the units of 
the square root of the string tension. The error given is statistical and 
does not include the systematic error. 



3.3 Masses of KK states 

Just like in the case of QCD4, we would like to analyze the behavior of the 
masses of the different KK modes. We will find very similar results: the KK 
modes wrapping the coordinate r are decoupling (even though a bit slower 
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10 



Figure 3.1: The dependence on a of the ratio r 



M, 



()++- 



M, 



0+H 



in QCD 3 . One 



can see that the ratio is very stable to changes in a, and reaches its 
asymptotic value quickly, a is given in units of u . 



than in QCD 4 ), while the other KK modes corresponding to states with 
angular momentum on the S 5 are not decoupling in the supergravity limit. 

First we consider the KK modes wrapping the compact r direction. Let 
us consider solutions to the Laplace equation V 2 $ = of the form 



$ = x(u] 



e ik-x e i/3T 



(3.12) 



The coordinate r is periodic with period Tjj 1 , where T# is given in eq. 
Therefore 

P = 2irT H n , (3.13) 



where n is an integer. Using the metric (|3.2j ) we find 



d u ((u 4 



Ur 



a 2 u 2 )ux'(u)) = u(-M 2 + n 2 4vr 2 T| 



u 



u 2 — a 2 r 

u 



) X (u). (3.14) 



One can see that just like in the case of QCD 4 there is an additional positive 
contribution to the masses, which grows like a 4 , therefore the masses of these 
KK states should grow as M\ K w a 4 . Thus these KK modes decouple from 
the spectrum, but slower than the corresponding KK modes in QCD4. The 



results of the numerical analysis are summarized in Table 3.2 and Figure 3.2 



Next we analyze the KK modes which correspond to states with angular 
momentum I = 1 on S 5 in the a = case. For a = these states have been 
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state 


value for a = 


value for a = 4 


KK 


5.79 


23.77 


KK* 


8.64 


24.63 


KK** 


11.50 


25.78 




14.36 


27.19 



Table 3.2: Masses of the KK modes wrapping the circle r in QCD 3 , 
using the same normalization as in Table p.l| . The first column gives the 
masses for a = while the second the masses for a = 4. Note that these 
states decouple quickly from the spectrum even in the supergravity 
approximation. 




Figure 3.2: The dependence on a of the ratio r = of the lowest 

++ glueball state compared to the KK states wrapping the r circle in 
units where uq = 1. These KK modes decouple from the spectrum in 
the supergravity approximation very quickly. 
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examined in [B], and found to be non-decoupling in the supergravity limit 
and including the lowest order a' corrections. Here we repeat this analysis 
and find (just like in the case of QCD 4 ) them to be non-decoupling even in 
the a — > oo case, in the supergravity limit. 

In order to do the analysis of these KK modes one needs to find the ex- 
plicit form of the spherical harmonics. The spherical harmonics on S d can 
be constructed in the following way. One takes S d embedded in R d+l , and 
expresses the Cartesian coordinates yi in terms of the angles. Then the spher- 
ical harmonics are just the functions C ilv " ,lk yi 1 ....yi k , where C is a symmetric 



traceless tensor p!J . This way, the simplest non-trivial spherical harmonic is 
just the coordinate yi itself. In the case of our QCD 3 theory, we actually have 
to use the "spheroidal coordinates" yi given in || on page 9. Thus one looks 
for solutions of the dilaton equation of the form f{u)e lk ' x y^ i = 1,2, 3, 4, 5, 6. 
For a = the isometry group of S 5 is 50(6), and the / = 1 KK mode is in 
the representation 6 of S0(6). Introducing the angular momentum a breaks 
SO (6) to SO (A) x U(l) x 17(1), and the 6 decomposes into 4+1 + 1. These 
states satisfy different eigenvalue equations. For i — 1, 2 (the two singlets 
are degenerate) the equation one gets is 

d u [u(u 4 ~ a 2 u 2 - - f(u) (k 2 u + 5u3 +jff~f \ = o. (3 . 15) 

Note that for a = this indeed reduces to the equation given in || for I = 1. 
We have numerically solved this equation, and find that the mass of these KK 
modes is growing slightly, until it becomes degenerate with the first excited 
state of the ++ glueball. Thus it does not decouple from the spectrum in 



the supergravity limit. The results are summarized Table [O] and Figure [373 
For the other 4 KK states which are in the 4 of SO (A) one finds the 
equation 

d u [u(u 4 - a 2 u 2 - l)f\u)} - f(u)(k 2 u + 5m 3 - 3a 2 u) = , (3.16) 

which for a = reproduces the equation in || with 1 = 1. The additional 
mass term is now negative, which actually makes these KK modes lighter in 
the a —>■ oo limit than for a = 0. However, they are still of the same order 
(and slightly heavier) than the ++ glueballs. The numerical results for this 



state are summarized in Table 3.4 and Figure 3.4 
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state 


value for a = 


value for a — > oo 


KK 


5.27 


7.05 


KK* 


8.29 


9.97 


KK** 


11.23 


12.87 




14.14 


15.76 



Table 3.3: Masses of the KK modes corresponding to the two degenerate 
singlet pieces of the I = 1 sextet of the original SO(Q) isometry in 
QCD 3 , using the same normalization as in Table |3J]. The first column 
gives the masses for a = while the second the masses in the a — > oo 
limit. Note that these states do not decouple from the spectrum in the 
supergravity approximation. 




Figure 3.3: The dependence of the ratio r = M K * + of the KK modes 
corresponding to the two singlet I = 1 states compared to the lowest ++ 
glueball state on the parameter a in units where uq — 1. This KK mode 
does not decouple from the spectrum in the supergravity approximation 
even in the a —>■ oo limit, even though it increases slightly and becomes 
exactly degenerate with the first excited glueball state ++ *. 
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state 


value for a = 


value for a — > oo 


KK 


5.27 


4.98 


KK* 


8.29 


8.14 


KK** 


11.23 


11.15 




14.14 


14.10 



Table 3.4: Masses of the KK modes corresponding to the quartet piece 
of the I = 1 sextet of the original SO (6) in QCD3, using the same 



normalization as in Table [3.11 The first column gives the masses for 
a = while the second the masses in the a — > 00 limit. Note that these 
states actually get lighter from a = to a = 00 in the supergravity 
approximation. 



1.3 



1.28 



1.26 



1.24 



1.22 



10 



of the KK states 

M o++ 



Figure 3.4: The dependence of the ratio r 
compared to the lowest ++ glueball state on the parameter a in units 
where Uq — 1. This KK mode does not decouple from the spectrum in 
the supergravity approximation even in the a — > 00 limit, they instead 
get even slightly lighter than for a = 0. 
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3.4 Free energy and Gluon condensation 

From the rotating D3 brane metric (see Eq. (3.16) in one can find the 
following formulas for the thermodynamic variables: 

m adm = ^m) 5 4 , nh2 m) = 7r 3 ; (317) 

47rG at 4 5 
S = 3 — — — 2mr#cosha, (3.18) 

Jh — — mi cosh a, (3.19) 

47tGat 



. /; _r H {2r 2 H + l 2 ) ^ _ (O0(j) 

47rmcosha ' 2m cosh a ' 



G, = ^ = 8<#r 6 (a') 4 , (3.21) 

2m cosh a sin a = Airg s Na' 2 . (3.22) 

One can check that they satisfy the first law of black hole thermodynamics 
Eq. (|2.31|) . In the limit a' — > (rescaling variables as in (|3.1| ) ), we get 

E = M ADM - M ext = ^V 3 N 2 u A , (3.23) 

T H = ^ , S = ±V 3 N 2 u H ul , (3.24) 
Tin Ztt 

Q H = i^f. J H = t -L Vs N 2 au 2 , (3.25) 

.2 2 U D ^ U 



where 



ir "-' r =4' B= u B w°-^ - (3 - 26) 



The free energy is then given by 



N 2 v 4 

E-T H S- n H J H = -V 3 —^ . (3.27) 



57T 2 



This gives for the gluon condensate the following expression: 

1 -Tr F 2 M) = -77^ = ^N 2 Bul . (3.28) 



¥m " w V 3 T H 87T 
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In terms of the Yang-Mills string tension, this is 



1 N 2 

"Tr FL(0)) = -— o 2 . (3.29) 



4#ym ^ " Stt A 

We find again that supergravity predicts that the gluon condensate is pro- 
portional to N 2 /\ times the string tension squared. The result expressed in 
terms of the string tension is thus independent of a. 



4 Conclusions 

In this paper we investigated quantitative aspects of large N SU (N) Yang- 
Mills theory in three and four dimensions using a one-parameter family of 
supergravity models related to non-extremal rotating D-branes. The new 
feature of this class of models is the decoupling of the KK modes associated 
with the compact D-brane coordinate as the angular momentum parameter 
is increased. The mass ratios for ordinary glueballs were found to be very 
stable against this variation. While the mass ratios of the ++ glueballs 
change only slightly compared to the case with zero angular momentum, 
there is a substantial change in the mass ratios of _+ , _+ * given in Eqs. 
( |2.20| ), ( |2.21| ), which for large a are in better agreement with the lattice values 



than for a = 0. 

It is worth emphasizing that the ratio a/uo should be large enough to 
have Mkk 3> ^giuebaii, but not infinite, since there are also string states 
winding around the compact D-brane coordinate with masses of order aRo 
that should decouple, i.e. M win( j ^> M g i ue b a ii. This requires Xu^/a 8 ^> 1, 
which is consistent with the condition that curvature invariants are small 
compared to the string scale ||. In general, for any given ratio a/uo which is 
large enough to decouple KK states from the low-energy physics, it is possible 
to choose A sufficiently large so that string winding states also decouple. 

We have found that the (5*0(3) or 5*0(4)) non-singlet KK modes with 
vanishing U(l) charge in the compact D-brane coordinate do not decouple in 
this class of models. One can hope that those KK modes may decouple in a 
model with more angular momenta (since there is room to take other limits). 
In this case the isometry group of the internal space is smaller. For example, 
in QCD 3 , for a = it is given by 50(6) x £7(1), whereas for a ^ it is 
50(4) x £7(1) x U(l). The isometry group of the model with the maximum 
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number of angular momenta only contains U(l) factors. This is consistent 
with the fact that in pure QCD there can only be singlets of the original 
R-symmetry. 

We have also found some features which seem to be universal, i.e. which 
do not depend on the extra supergravity parameter. In particular, both in 
QCD 3 and QCD 4 supergravity gives a gluon condensate of the form 3-0" 2 , 
with a coefficient which is the same for all models parametrized by a. Another 
feature that seems to be common to all supergravity models is a topological 
susceptibility of the form Act 2 , with a coefficient which is independent of N 
but depends on a/uo- This result suggests that in the regime A ^> 1 the 
r)' particle of QCD 4 with N = 3 is much heavier than other mesons (whose 
masses are proportional to the string tension). 

Acknowledgements 

We thank M. Dine, A. Hashimoto, K. Hori, K. Jansen, C. Morningstar, H. 
Ooguri, M. Peardon, and A. Zaffaroni for useful discussions. C.C. is a re- 
search fellow of the Miller Institute for Basic Research in Science. C.C. and 
J.T. are supported in part the U.S. Department of Energy under Contract 
DE-AC03-76SF00098, and in part by the National Science Foundation un- 
der grant PHY-95- 14797. The work of J.R. is supported by the European 
Commission TMR programme grant ERBFMBI-CT96-0982. 

References 

[1] J. M. Maldacena, "The Large N Limit of Superconformal Field Theories 
and Supergravity", hep-th/9711200. 

[2] E. Witten, "Anti-de Sitter space, thermal phase transition, and confine- 
ment in gauge theories", hep-th/9803131. 

[3] D. J. Gross and H. Ooguri, "Aspects of large N gauge theory dynamics as 
seen by string theory" , Phys. Rev. D58 (1998) 106002, hep-th/9805129. 

[4] C. Csaki, H. Ooguri, Y. Oz, and J. Terning, "Glueball mass spectrum 
from supergravity", hep-th/9806021. 



28 



[5] R. de Mello Koch, A. Jevicki, M. Mihailescu and J. Nunes, "Evaluation 
of glueball masses from supergravity" , Phys. Rev. D58 (1998) 105009, 
hep-th/9806125; M. Zyskin, "A note on the glueball mass spectrum", 
hep-th/9806128. 

[6] H. Ooguri, H. Robins and J. Tannenhauser, "Glueballs and their Kaluza- 
Klein cousins", hep-th/9806171. 

[7] S. S. Gubser, I. R. Klebanov, and A. M. Polyakov, "Gauge theory cor- 
relators from noncritical string theory," Phys. Lett. B428 (1998) 105, 
hep-th/9802109. 

[8] E. Witten, "Anti-de Sitter space and holography," hep-th/9802150. 

[9] J.G. Russo, "New compactifications of supergravities and large N 
QCD", hep-th/9808117. 

[10] M. J. Teper, "Physics from the lattice: Glueballs in QCD: Topology: 
SU(N) for all N," hep-lat/9711011. 

[11] C. Morningstar and M. Peardon, "Efficient glueball simulations on 
anisotropic lattices," Phys. Rev. D56 (1997) 4043, hep-lat/9704011. 

[12] A. Hashimoto and Y. Oz, "Aspects of QCD dynamics from string the- 
ory", hep-th/9809106. 

[13] M. Peardon, "Coarse lattice results for glueballs and hybrids," 
Nucl. Phys. B (Proc. Suppl.) 63, 22 (1998), hep-lat/9710029; C. Morn- 
ingstar and M. Peardon, in preparation. 

[14] C. Michael, "Hadronic spectroscopy from the lattice: Glueballs and Hy- 
brid Mesons", hep-ph/9810415. 

[15] S. Dalley and B. van de Sande, "Glueball calculations in large N c gauge 
theory", hep-ph/9810236. 

[16] G. Horowitz and A. Sen, "Rotating black holes which saturate a Bogo- 
mol'nyi bound", Phys. Rev. D53 (1996) 808, hep-th/9509108. 

[17] M. Cvetic and D. Youm, "Rotating intersecting M-branes", Nucl. Phys. 
B499 (1997) 253, hep-th/9612229. 



29 



[18] I. Klebanov and A. A Tseytlin, "Entropy of near-extremal black p- 
branes", Nucl. Phys. B475 (1996) 164, hep-th/9604089. 

[19] E. Witten, "Current algebra theorems for the U(l) Goldstone boson," 
Nucl. Phys. B156 (1979) 269. 

[20] G. Veneziano, U U(1) without instantons," Nucl. Phys. B159 (1979) 213. 

[21] E. Witten, "Theta dependence in the large N limit of four- dimensional 
gauge theories," Phys. Rev. Lett. 81 (1998) 2862, hep-th/9807109. 

[22] J.M. Maldacena, "Wilson loops in large N field theories", 
Phys. Rev. Lett. 80 (1998) 4859, hep-th/9803002. 

[23] A. Brandhuber, N. Itzhaki, J. Sonnenschein and S. Yankielowicz, 
"Wilson loops, confinement and phase transitions in large N gauge 
theories from supergravity" , J. High Energy Phys. 06 (1998) 001, 
hep-th/9803263. 

[24] S. Rey and J. Yee, "Macroscopic Strings as Heavy Quarks in Large 
N Gauge Theory and Anti-de Sitter Supergravity", hep-th/9803001; 
S. Rey, S. Theisen and J. Yee, "Wilson-Polyakov Loop at Finite Tem- 
perature in Large N Gauge Theory and Anti-de Sitter Supergravity", 
Nucl. Phys. B527 (1998) 171, hep-th/9803135. 

[25] M. J. Teper, "SU(N) gauge theories in (2 + 1) dimensions", 
hep-lat/9804008. 

[26] S. Lee, S. Minwalla, M. Rangamani and N. Seiberg, "Three point 
functions of chiral operators in D = 4, N = 4 SYM at large JV", 
hep-th/9806074. 



30 



